Understanding conceptus-maternal interactions: what tools do we need to develop?

Communication between the maternal endometrium and developing embryo/conceptus is critical to support successful pregnancy to term. Studying the peri-implantation period of pregnancy is critical as this is when most pregnancy loss occurs in cattle. Our current understanding of these interactions is limited, due to the lack of appropriate in vitro models to assess these interactions. The endometrium is a complex and heterogeneous tissue that is regulated in a transcriptional and translational manner throughout the oestrous cycle. While there are in vitro models to study endometrial function, they are static and 2D in nature or explant models and are limited in how well they recapitulate the in vivo endometrium. Recent developments in organoid systems, microfluidic approaches, extracellular matrix biology, and in silico approaches provide a new opportunity to develop in vitro systems that better model the in vivo scenario. This will allow us to investigate in a more high-throughput manner the fundamental molecular interactions that are required for successful pregnancy in cattle.

Caecilian Genomes Reveal the Molecular Basis of Adaptation and Convergent Evolution of Limblessness in Snakes and Caecilians.

We present genome sequences for the caecilians Geotrypetes seraphini (3.8 Gb) and Microcaecilia unicolor (4.7 Gb), representatives of a limbless, mostly soil-dwelling amphibian clade with reduced eyes, and unique putatively chemosensory tentacles. More than 69% of both genomes are composed of repeats, with retrotransposons being the most abundant. We identify 1,150 orthogroups that are unique to caecilians and enriched for functions in olfaction and detection of chemical signals. There are 379 orthogroups with signatures of positive selection on caecilian lineages with roles in organ development and morphogenesis, sensory perception, and immunity amongst others. We discover that caecilian genomes are missing the zone of polarizing activity regulatorysequence (ZRS) enhancer of Sonic Hedgehog which is also mutated in snakes. In vivo deletions have shown ZRS is required for limb development in mice, thus, revealing a shared molecular target implicated in the independent evolution of limblessness in snakes and caecilians.

A burst of genomic innovation at the origin of placental mammals mediated embryo implantation.

The origin of embryo implantation in mammals ~148 million years ago was a dramatic shift in reproductive strategy, yet the molecular changes that established mammal implantation are largely unknown. Although progesterone receptor signalling predates the origin of mammals and is highly conserved in, and critical for, successful mammal pregnancy, it alone cannot explain the origin and subsequent diversity of implantation strategies throughout the placental mammal radiation. MiRNAs are known to be flexible and dynamic regulators with a well-established role in the pathophysiology of mammal placenta. We propose that a dynamic core microRNA (miRNA) network originated early in placental mammal evolution, responds to conserved mammal pregnancy cues (e.g. progesterone), and facilitates species-specific responses. Here we identify 13 miRNA gene families that arose at the origin of placental mammals and were subsequently retained in all descendent lineages. The expression of these miRNAs in response to early pregnancy molecules is regulated in a species-specific manner in endometrial epithelia of species with extreme implantation strategies (i.e. bovine and human). Furthermore, this set of miRNAs preferentially target proteins under positive selective pressure on the ancestral eutherian lineage. Discovery of this core embryo implantation toolkit and specifically adapted proteins helps explain the origin and evolution of implantation in mammals.

MicroRNAs emerging coordinate with placental mammals alter pathways in endometrial epithelia important for endometrial function.

We tested the hypothesis that conserved placental mammal-specific microRNAs and their targets facilitate endometrial receptivity to implantation. Expression of miR-340-5p, -542-3p, and -671-5p was regulated by exposure of endometrial epithelial cells to progesterone (10 µg/ml) for 24 h coordinate with 1,713 of their predicted targets. Proteomic analysis of cells transfected with miRNA mimic/inhibitor (48 h: n = 3) revealed 1,745 proteins altered by miR-340-5p (mimic; 1,369, inhibitor; 376) of which 171 were predicted targets and P4-regulated. MiR-542-3p altered 2,353 (mimic; 1,378, inhibitor; 975) 100 which were mirDB predicted, including 46 P4-regulated. MiR-671-5p altered 1,744 proteins (mimic; 1,252, inhibitor; 492) 95 of which were predicted targets and 46 P4-regulated. All miRNAs were detected in luteal phase endometrial biopsies, irrespective of pregnancy outcomes. miR-340-5p expression increased in biopsies from individuals suffering previous and subsequent miscarriage compared to those with subsequent live birth. Dysfunction of these miRNAs and their targets contribute to endometrial-derived recurrent pregnancy loss.

Improving Orthologous Signal and Model Fit in Datasets Addressing the Root of the Animal Phylogeny.

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference.

Filtering artifactual signal increases support for Xenacoelomorpha and Ambulacraria sister relationship in the animal tree of life.

Conflicting studies place a group of bilaterian invertebrates containing xenoturbellids and acoelomorphs, the Xenacoelomorpha, as either the primary emerging bilaterian phylum1,2,3,4,5,6 or within Deuterostomia, sister to Ambulacraria.7,8,9,10,11 Although their placement as sister to the rest of Bilateria supports relatively simple morphology in the ancestral bilaterian, their alternative placement within Deuterostomia suggests a morphologically complex ancestral bilaterian along with extensive loss of major phenotypic traits in the Xenacoelomorpha. Recent studies have questioned whether Deuterostomia should be considered monophyletic at all.10,12,13 Hidden paralogy and poor phylogenetic signal present a major challenge for reconstructing species phylogenies.14,15,16,17,18 Here, we assess whether these issues have contributed to the conflict over the placement of Xenacoelomorpha. We reanalyzed published datasets, enriching for orthogroups whose gene trees support well-resolved clans elsewhere in the animal tree.16 We find that most genes in previously published datasets violate incontestable clans, suggesting that hidden paralogy and low phylogenetic signal affect the ability to reconstruct branching patterns at deep nodes in the animal tree. We demonstrate that removing orthogroups that cannot recapitulate incontestable relationships alters the final topology that is inferred, while simultaneously improving the fit of the model to the data. We discover increased, but ultimately not conclusive, support for the existence of Xenambulacraria in our set of filtered orthogroups. At a time when we are progressing toward sequencing all life on the planet, we argue that long-standing contentious issues in the tree of life will be resolved using smaller amounts of better quality data that can be modeled adequately.19.

Ribosome heterogeneity in Drosophila melanogaster gonads through paralog-switching.

Ribosomes have long been thought of as homogeneous macromolecular machines, but recent evidence suggests they are heterogeneous and could be specialised to regulate translation. Here, we have characterised ribosomal protein heterogeneity across 4 tissues of Drosophila melanogaster. We find that testes and ovaries contain the most heterogeneous ribosome populations, which occurs through a combination of paralog-enrichment and paralog-switching. We have solved structures of ribosomes purified from in vivo tissues by cryo-EM, revealing differences in precise ribosomal arrangement for testis and ovary 80S ribosomes. Differences in the amino acid composition of paralog pairs and their localisation on the ribosome exterior indicate paralog-switching could alter the ribosome surface, enabling different proteins to regulate translation. One testis-specific paralog-switching pair is also found in humans, suggesting this is a conserved site of ribosome heterogeneity. Overall, this work allows us to propose that mRNA translation might be regulated in the gonads through ribosome heterogeneity, providing a potential means of ribosome specialisation.

Adaptive Evolution in TRIF Leads to Discordance between Human and Mouse Innate Immune Signaling.

The TIR domain-containing adapter inducing IFN-ß (TRIF) protein is an innate immune system protein that mediates the MyD88-independent toll-like receptor response pathway in mice and humans. Previously, we identified positive selection at seven distinct residues in mouse TRIF (mTRIF), as compared with human and other mammalian orthologs, thus predicting protein functional shift in mTRIF. We reconstructed TRIF for the most recent common ancestor of mouse and human, and mutated this at the seven sites to their extant mouse/human states. We overexpressed these TRIF mutants in immortalized human and mouse cell lines and monitored TRIF-dependent cytokine production and gene expression induction. We show that optimal TRIF function in human and mouse is dependent on the identity of the positively selected sites. These data provide us with molecular data relating observed differences in response between mouse and human MyD88-independent signaling in the innate immune system with protein functional change.

Cytoplasmic long noncoding RNAs are differentially regulated and translated during human neuronal differentiation.

The expression of long noncoding RNAs is highly enriched in the human nervous system. However, the function of neuronal lncRNAs in the cytoplasm and their potential translation remains poorly understood. Here we performed Poly-Ribo-Seq to understand the interaction of lncRNAs with the translation machinery and the functional consequences during neuronal differentiation of human SH-SY5Y cells. We discovered 237 cytoplasmic lncRNAs up-regulated during early neuronal differentiation, 58%-70% of which are associated with polysome translation complexes. Among these polysome-associated lncRNAs, we find 45 small ORFs to be actively translated, 17 specifically upon differentiation. Fifteen of 45 of the translated lncRNA-smORFs exhibit sequence conservation within Hominidea, suggesting they are under strong selective constraint in this clade. The profiling of publicly available data sets revealed that 8/45 of the translated lncRNAs are dynamically expressed during human brain development, and 22/45 are associated with cancers of the central nervous system. One translated lncRNA we discovered is LINC01116, which is induced upon differentiation and contains an 87 codon smORF exhibiting increased ribosome profiling signal upon differentiation. The resulting LINC01116 peptide localizes to neurites. Knockdown of LINC01116 results in a significant reduction of neurite length in differentiated cells, indicating it contributes to neuronal differentiation. Our findings indicate cytoplasmic lncRNAs interact with translation complexes, are a noncanonical source of novel peptides, and contribute to neuronal function and disease. Specifically, we demonstrate a novel functional role for LINC01116 during human neuronal differentiation.

Spectrum of pathogenic variants and founder effects in amelogenesis imperfecta associated with MMP20.

Amelogenesis imperfecta (AI) describes a heterogeneous group of developmental enamel defects that typically have Mendelian inheritance. Exome sequencing of 10 families with recessive hypomaturation AI revealed four novel and one known variants in the matrix metallopeptidase 20 (MMP20) gene that were predicted to be pathogenic. MMP20 encodes a protease that cleaves the developing extracellular enamel matrix and is necessary for normal enamel crystal growth during amelogenesis. New homozygous missense changes were shared between four families of Pakistani heritage (c.625G>C; p.(Glu209Gln)) and two of Omani origin (c.710C>A; p.(Ser237Tyr)). In two families of UK origin and one from Costa Rica, affected individuals were homozygous for the previously reported c.954-2A>T; p.(Ile319Phefs*19) variant. For each of these variants, microsatellite haplotypes appeared to exclude a recent founder effect, but elements of haplotype were conserved, suggesting more distant founding ancestors. New compound heterozygous changes were identified in one family of the European heritage: c.809_811+12delinsCCAG; p.(?) and c.1122A>C; p.(Gln374His). This report further elucidates the mutation spectrum of MMP20 and the probable impact on protein function, confirms a consistent hypomaturation phenotype and shows that mutations in MMP20 are a common cause of autosomal recessive AI in some communities.

The role of CAPG in molecular communication between the embryo and the uterine endometrium: Is its function conserved in species with different implantation strategies?

During the preimplantation period of pregnancy in eutherian mammals, transcriptional and proteomic changes in the uterine endometrium are required to facilitate receptivity to an implanting blastocyst. These changes are mediated, in part, by proteins produced by the developing conceptus (inner cell mass and extraembryonic membranes). We hypothesized that this common process in early pregnancy in eutheria may be facilitated by highly conserved conceptus-derived proteins such as macrophage capping protein (CAPG). We propose that CAPG may share functionality in modifying the transcriptome of the endometrial epithelial cells to facilitate receptivity to implantation in species with different implantation strategies. A recombinant bovine form of CAPG (91% sequence identity between bovine and human) was produced and bovine endometrial epithelial (bEECs) and stromal (bESCs) and human endometrial epithelial cells (hEECs) were cultured for 24 hours with and without recombinant bovine CAPG (rbCAPG). RNA sequencing and quantitative real-time PCR analysis were used to assess the transcriptional response to rbCAPG (Control, vehicle, CAPG 10, 100, 1000 ng/mL: n = 3 biological replicates per treatment per species). Treatment of bEECs with CAPG resulted in alterations in the abundance of 1052 transcripts (629 increased and 423 decreased) compared to vehicle controls. Treatment of hEECs with bovine CAPG increased expression of transcripts previously known to interact with CAPG in different systems (CAPZB, CAPZA2, ADD1, and ADK) compared with vehicle controls (P < .05). In conclusion, we have demonstrated that CAPG, a highly conserved protein in eutherian mammals, elicits a transcriptional response in the endometrial epithelium in species with different implantation strategies that may contribute to pregnancy success.

Protection against SIV in Rhesus Macaques Using Albumin and CD4-Based Vector-Mediated Gene Transfer.

Antibody-like molecules were evaluated with potent simian immunodeficiency virus (SIV) neutralizing properties (immunoadhesins) that were delivered by a recombinant adeno-associated virus (rAAV) vector in the SIV-infected rhesus macaque model. When injected intramuscularly into the host, the vector directs in vivo production of the transgenes with antibody-like binding properties that lead to serum neutralizing activity against SIV. To extend the half-life of the immunoadhesins, rhesus cluster of differentiation 4 (CD4) and a single-chain antibody (4L6) were fused with albumin molecules, and these constructs were tested in our model of SIV infection. Antibody-based immunoadhesins provided high serum neutralizing titers against the original SIV strain. CD4-based immunoadhesins provided a wider spectrum of neutralization against different SIV strains in comparison to antibody-based therapeutics and had the potential to protect against high viral challenging doses. Although the albumin-antibody fusion immunoadhesin provided strong and prolonged protection of the immunized animals against SIV challenge, the albumin-CD4 fusion altered the specificity and decreased the overall protection effectiveness of the immunoadhesin in comparison to the antibody-based molecules. Albumin-based immunoadhesins increase in vivo longevity of the immune protection; however, they present challenges likely linked to the induction of anti-immunoadhesin antibodies.

Atlanto-Occipital Transarticular Screw Fixation for the Treatment of Traumatic Occipitocervical Instability in the Pediatric Population.

BACKGROUND: Atlanto-occipital transarticular screw fixation (AOTSF) has rarely been reported for fixation of the craniovertebral junction (CVJ). METHODS: A retrospective chart review of all pediatric patients (less than 18 years of age) with an attempt of AOTSF for fixation of traumatic CVJ instability was conducted. RESULTS: A total of 4 patients (2 boys and 2 girls; ages 2, 3, 5, and 8 years) who suffered from acute traumatic CVJ instability managed during 2007-2018 underwent an attempted AOTSF. In 2 patients, this method was technically successful. In the other 2 instances, we were not able to engage the screw into the occipital condyle. These were converted to standard occipital plate, rod, and screw fixation. All were placed in a halo subsequently for a minimum of 3 months. Three patients were fused at last follow-up (range, 17-48 months). One patient after successful AOTSF did not fuse. There were no surgical complications or revision procedures. CONCLUSIONS: AOTSF was feasible in half of pediatric patients suffering from traumatic CVJ instability. Therefore, intraoperative salvage options and strategies should be on hand readily. In the pediatric population, where bony anatomy may pose challenges to fixation, this technique may offer a viable first-line option in selected cases, despite the overall modest success rate.

New missense variants in RELT causing hypomineralised amelogenesis imperfecta.

Amelogenesis imperfecta (AI) is a heterogeneous group of genetic diseases characterised by dental enamel malformation. Pathogenic variants in at least 33 genes cause syndromic or non-syndromic AI. Recently variants in RELT, encoding an orphan receptor in the tumour necrosis factor (TNF) superfamily, were found to cause recessive AI, as part of a syndrome encompassing small stature and severe childhood infections. Here we describe four additional families with autosomal recessive hypomineralised AI due to previously unreported homozygous mutations in RELT. Three families carried a homozygous missense variant in the fourth exon (c.164C>T, p.(T55I)) and a fourth family carried a homozygous missense variant in the 11th exon (c.1264C>T, p.(R422W)). We found no evidence of additional syndromic symptoms in affected individuals. Analyses of tooth microstructure with computerised tomography and scanning electron microscopy suggest a role for RELT in ameloblasts' coordination and interaction with the enamel matrix. Microsatellite genotyping in families segregating the T55I variant reveals a shared founder haplotype. These findings extend the RELT pathogenic variant spectrum, reveal a founder mutation in the UK Pakistani population and provide detailed analysis of human teeth affected by this hypomineralised phenotype, but do not support a possible syndromic presentation in all those with RELT-variant associated AI.

Strigolactone synthesis is ancestral in land plants, but canonical strigolactone signalling is a flowering plant innovation.

BACKGROUND: Strigolactones (SLs) are an important class of carotenoid-derived signalling molecule in plants, which function both as exogenous signals in the rhizosphere and as endogenous plant hormones. In flowering plants, SLs are synthesized by a core pathway of four enzymes and are perceived by the DWARF14 (D14) receptor, leading to degradation of SMAX1-LIKE7 (SMXL7) target proteins in a manner dependent on the SCFMAX2 ubiquitin ligase. The evolutionary history of SLs is poorly understood, and it is not clear whether SL synthesis and signalling are present in all land plant lineages, nor when these traits evolved. RESULTS: We have utilized recently-generated genomic and transcriptomic sequences from across the land plant clade to resolve the origin of each known component of SL synthesis and signalling. We show that all enzymes in the core SL synthesis pathway originated at or before the base of land plants, consistent with the previously observed distribution of SLs themselves in land plant lineages. We also show that the late-acting enzyme LATERAL BRANCHING OXIDOREDUCTASE (LBO) may be considerably more ancient than previously thought. We perform a detailed phylogenetic analysis of SMXL proteins and show that specific SL target proteins only arose in flowering plants. We also assess diversity and protein structure in the SMXL family, identifying several previously unknown clades. CONCLUSIONS: Overall, our results suggest that SL synthesis is much more ancient than canonical SL signalling, consistent with the idea that SLs first evolved as rhizosphere signals and were only recruited much later as hormonal signals.

Gene Fusions Derived by Transcriptional Readthrough are Driven by Segmental Duplication in Human.

Gene fusion occurs when two or more individual genes with independent open reading frames becoming juxtaposed under the same open reading frame creating a new fused gene. A small number of gene fusions described in detail have been associated with novel functions, for example, the hominid-specific PIPSL gene, TNFSF12, and the TWE-PRIL gene family. We use Sequence Similarity Networks and species level comparisons of great ape genomes to identify 45 new genes that have emerged by transcriptional readthrough, that is, transcription-derived gene fusion. For 35 of these putative gene fusions, we have been able to assess available RNAseq data to determine whether there are reads that map to each breakpoint. A total of 29 of the putative gene fusions had annotated transcripts (9/29 of which are human-specific). We carried out RT-qPCR in a range of human tissues (placenta, lung, liver, brain, and testes) and found that 23 of the putative gene fusion events were expressed in at least one tissue. Examining the available ribosome foot-printing data, we find evidence for translation of three of the fused genes in human. Finally, we find enrichment for transcription-derived gene fusions in regions of known segmental duplication in human. Together, our results implicate chromosomal structural variation brought about by segmental duplication with the emergence of novel transcripts and translated protein products.

Return to the Sea, Get Huge, Beat Cancer: An Analysis of Cetacean Genomes Including an Assembly for the Humpback Whale (Megaptera novaeangliae).

Cetaceans are a clade of highly specialized aquatic mammals that include the largest animals that have ever lived. The largest whales can have ∼1,000× more cells than a human, with long lifespans, leaving them theoretically susceptible to cancer. However, large-bodied and long-lived animals do not suffer higher risks of cancer mortality than humans-an observation known as Peto's Paradox. To investigate the genomic bases of gigantism and other cetacean adaptations, we generated a de novo genome assembly for the humpback whale (Megaptera novaeangliae) and incorporated the genomes of ten cetacean species in a comparative analysis. We found further evidence that rorquals (family Balaenopteridae) radiated during the Miocene or earlier, and inferred that perturbations in abundance and/or the interocean connectivity of North Atlantic humpback whale populations likely occurred throughout the Pleistocene. Our comparative genomic results suggest that the evolution of cetacean gigantism was accompanied by strong selection on pathways that are directly linked to cancer. Large segmental duplications in whale genomes contained genes controlling the apoptotic pathway, and genes inferred to be under accelerated evolution and positive selection in cetaceans were enriched for biological processes such as cell cycle checkpoint, cell signaling, and proliferation. We also inferred positive selection on genes controlling the mammalian appendicular and cranial skeletal elements in the cetacean lineage, which are relevant to extensive anatomical changes during cetacean evolution. Genomic analyses shed light on the molecular mechanisms underlying cetacean traits, including gigantism, and will contribute to the development of future targets for human cancer therapies.

Susceptibility to SIV Infection After Adenoviral Vaccination in a Low Dose Rhesus Macaque Challenge Model.

BACKGROUND: Vaccination with the Merck human adenovirus serotype-5 (HAdV-5) vectored HIV-1 subtype B gag/pol/nef vaccine was unexpectedly associated with enhanced susceptibility to HIV-1 infection in uncircumcised HAdV-5 seropositive men. It has been hypothesized that vaccination may have resulted in activated CD4+ T lymphocytes trafficking to mucosal sites thereby increasing targets for HIV infection. We have previously shown that AdV-vector vacci-nation in rhesus macaques resulted in an increase in the frequency of activated mucosal CD4+ T cells. However, whether this increase in activation is sufficient to increase susceptibility to HIV/SIV infection is unclear. METHODS: To examine this scenario, we developed a preliminary, proof-of-concept vaccination-challenge model in order to examine vaccine-induced SIV susceptibility in rhesus macaques. Rhesus macaques (n = 10/group) were vaccinated with a simian AdV-7 (SAdV-7)-vector encoding an irrelevant insert (SARS spike) and challenged 5 weeks post-prime in an escalating dosing regimen starting with sub-infectious doses (1:10,000 or 2TCID50) of SIVmac251. RESULTS: In contrast to our previous study, the SAdV-7 vaccine regimen did not induce detectable mucosal CD4+ T cell activation at the time points assessed in animals obtained from a different vendor and housed in a different facility. Within the power of the study, we did not observe significantly increased SIV acquisition in SAdV-7-vaccinated (5/10) versus placebo-vaccinated (3/10) macaques after repeated low-dose intra-rectal SIVmac251 challenge (P < 0.2). CONCLUSIONS: These results lay groundwork for future experiments to assess vaccine-induced SIV susceptibility in rhesus macaques. Further larger-scale studies are necessary to confirm the AdV-vector vaccination associated trend towards increased SIV/HIV acquisition and clarify associated mechanisms.

Paternally Expressed Imprinted Genes under Positive Darwinian Selection in Arabidopsis thaliana.

Genomic imprinting is an epigenetic phenomenon where autosomal genes display uniparental expression depending on whether they are maternally or paternally inherited. Genomic imprinting can arise from parental conflicts over resource allocation to the offspring, which could drive imprinted loci to evolve by positive selection. We investigate whether positive selection is associated with genomic imprinting in the inbreeding species Arabidopsis thaliana. Our analysis of 140 genes regulated by genomic imprinting in the A. thaliana seed endosperm demonstrates they are evolving more rapidly than expected. To investigate whether positive selection drives this evolutionary acceleration, we identified orthologs of each imprinted gene across 34 plant species and elucidated their evolutionary trajectories. Increased positive selection was sought by comparing its incidence among imprinted genes with nonimprinted controls. Strikingly, we find a statistically significant enrichment of imprinted paternally expressed genes (iPEGs) evolving under positive selection, 50.6% of the total, but no such enrichment for positive selection among imprinted maternally expressed genes (iMEGs). This suggests that maternally- and paternally expressed imprinted genes are subject to different selective pressures. Almost all positively selected amino acids were fixed across 80 sequenced A. thaliana accessions, suggestive of selective sweeps in the A. thaliana lineage. The imprinted genes under positive selection are involved in processes important for seed development including auxin biosynthesis and epigenetic regulation. Our findings support a genomic imprinting model for plants where positive selection can affect paternally expressed genes due to continued conflict with maternal sporophyte tissues, even when parental conflict is reduced in predominantly inbreeding species.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics.

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.